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ABSTRACT 

Observations of the ELAIS-Nl field taken at 610 MHz with the Giant Metrewave Radio Tele- 
scope are presented. Nineteen pointings were observed, covering a total area of ^ 9 deg^ with 
a resolution of 6 x 5 arcsec^, PA +45°. Four of the pointings were deep observations with an 
rms of ^ 40 /X Jy before primary beam correction, with the remaining fifteen pointings having 
an rms of ^ 70 /xJy. The techniques used for data reduction and production of a mosaicked 
image of the region are described, and the final mosaic is presented, along with a catalogue 
of 2500 sources detected above 6(7. This work complements the large amount of optical and 
infrared data already available on the region. We calculate 610-MHz source counts down to 
270 ixiy, and find further evidence for the turnover in differential number counts below 1 mJy, 
previously seen at both 610 MHz and 1.4 GHz. 
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1 INTRODUCTION 

The Spitzer Wide -area Infrared Extragalactic (SWIRE; 
iLonsdale et id] l2003h survey has the largest sky coverage of 
the legacy surveys b eing performed by the Spitzer Space Telescope 
( IWemeretal1l2004h . A total area of '~49 deg^ of sky has been 
observed with the Infrared Array Camera (IRAC; Fazio et al. 2004) 
and M ultiband Imaging Photometer for Spitzer (MIPS; Reike et al.| 
l2004h instruments at 3.6, 4.5, 5.8, 8, 24, 70 and 160 ^m. The 
survey is broken down into six tields, three in the northern sky - 
ELAIS-Nl, ELAIS-N2 and the Lockman Hole - and three in the 
south - ELAIS-S 1 , Chandra Deep Field South and the XMM-Large 
Scale Structure (XMM-LSS) field. All six regions were selected to 
be away from the Galactic disk, in order to minimize background 
cirrus emission. 

There is a large amount of multi-wavelength information 
available on all six SWIRE fields. The three ELAIS fie lds were ob- 
serve d as part of the European Large- Area ISO Survey jOliver et al.l 
which also included another northern (-N3) and southern (- 
S2) field. The Infrared Space Observatory (ISO) observed these 
regions at 6.7, 15, 90 and 175 /im, and a large number of fol- 
lowup observations were carried out in the optical, infrared and 
radio bands. A band-merged catalogue, containing the ISO data, 
along with U, g' , r', i', Z, y, H and AT-band detections, and radio 
observations at 1.4 GHz has been produced - for more details, see 
[Powan-Robinson et al. (2004), and references therein. 

Observations of the ELAIS-Nl region were taken with 
Spitzer in 2004 January, covering ~ 9 deg^ with the IRAC and 
MIPS instruments. The source catalogues have been produced, 
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and are available online jSurace et al. I l2004l) . containing over 
280,000 sources. The UK Infrared Deep Sky Survey (UKIDSS; 
[Lawrence et al.l l2007h intends to cover the ELAIS-Nl region in 
its Deep Extragalactic Survey plan, observing the full field in 
the J, H and /f-bands to a depth of A" = 21 mag. This will be 
a great improvement over the currently available surveys, which 
have a sensitivity li mit of ~ 1 8 mag in the K band. Data Release 2 
dWarren et al.l2007t) of UKIDSS contains early shallow data on the 
ELAI S-Nl region. Furt her surve ys have been carried out in the R- 
band l lFadda et al]|2004l) . i n Hg (Pascual et al. 2001), and with the 
Chandra X-ray telescope ( [Manners e t al. 2003; Fra nceshini et al.l 
There have been several red shift surveys of the region 
dTrichas et alj|2006l ; Iserta et al.|[2007l) . and the ELAIS-Nl region 
was also partially covered by t he Sloan Digital Sky Survey (SDSS; 
lAdelman-McCartv et alj|2007l) . 

While there have been a great number of observations of 
the ELAIS-Nl region at optical and infrared wavelengths, there is 
comparatively little radio information available. The existing VLA 
1.4 G Hz survey of the three northern ELAIS fields jCiliegi et al.l 
11999"), which has been included into the band-merged catalogue 
of Rowan-Robinson et al. (2004), reaches a 5(T limit of 0.135 mJy 
over 0.12 deg^ but only a 1.15 mJy limit over its f ull coverage 
area of 4.22 deg^ The NVSS dCondon et alJll998t) and FIRST 
( lBeckeretal1ll995h surveys both cover the ELAIS-Nl region, but 
only to relatively shallow 5a limits of 2.25 and 0.75 mJy respec- 
tively. A recent study of polarised compact sources (Taylor et al] 
l2007h at 1420 MHz is underway, using the Dominion Radio As- 
trophysical Observatory Synthesis Telescope (DRAO ST) cen- 
tered on le'^ll^'OO', +55°00'00" and covering 7.4 deg^. The first 
30 per cent of observations have been completed, with maps in 
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Stokes I, Q and U being produced with a maximum sensitivity of 
78 /xJy beam^', although with a resolution of ~1 arcmin^. 

In order to extend the information on this region, a much larger 
deep radio survey is required. In this paper, we present observations 
of the ELAIS-Nl survey field taken at 610 MHz with the Giant Me- 
trewave Radio Telescope (GMRT; Ananthakrishnan 2005i) , cover- 
ing ~ 9 deg- of sky with a resolution of 6 x 5 arcsec^, PA +45°, 
centred on 16^11™00% +55°00'00" (J2000 coordinates, which are 
used throughout this paper). This survey, in combination with the 
deep Spitzer data, will be used to s tudy the infrare d/radio correla- 
tion for star-forming systems (e.g. lAppleton et al] |2004). and the 
link between the triggering of star formation and AGN activity, as 
well as the properties of the faint radio population at 610 MHz. 

In Section |2] we describe the observations and data reduction 
techniques used in the creation of the survey. Section [3] presents 
the mosaic and a source catalogue containing 2500 sources above 
6<T, along with a sample of extended sources. In Section|4]we con- 
struct the 610 MHz differential source counts, and compare them 
to previous works. 



2 OBSERVATIONS AND DATA REDUCTION 

The ELAIS-Nl region was observed for 25 hours, spread over three 
days in 2003 August, using the GMRT operating at 610 MHz. Nine- 
teen pointings were observed, centred on 16''H™00% +55°00'00" 
and spaced by 36' in a hexagonal grid (as shown in Fig.[Tll in order 
to get nearly uniform coverage over the region. Two 16 MHz side- 
bands were observed, each split into 128 spectral channels, with a 
16.9 s integration time. 

The flux density scale was set through observations of 3C48 
or 3C286, at the beginning and end of each observing session. The 
MPS task SETJY was used to calculate 610 MHz flux densities of 
29. 4 and 21.1 Jy, res pectively, using the AIPS implementation of 
the lSaars et alj ( Il977h scale. Each field was observed for a series of 
interleaved 9 min scans in order to maximise the uv coverage, and a 
nearby phase calibrator, J1634+627, was observed for four minutes 
between every three scans to monitor any time-dependent phase 



and amplitude fluctuations of the telescope. The measured phase 
typically varied by less than 10 degrees between phase calibrator 
observations. Most of the fields were observed for four scans (36 
mins in total), while four fields, B4, C3, C5 and D4, were observed 
for significantly longer, in order to give a deeper region within the 
main survey as there is deep optical inform ation available on this 
region from UKIDSS jWarren et ai]|2007h . B4 and D4 were ob- 
served for 1 1 scans (99 mins) while C3 and C5 were observed for 
12 scans (108 mins). 

The imaging strategy was very similar to that used for our 
610 MHz GMRT observ ations of the Sp itzer extragalactic First 
Look Survey field (xFLS: lGam et alj|2007h . An error in the times- 
tamps o f the uv data wa s corrected using UVFXT (for further de- 
tails see Gam et al .120071) . This error is only found for observations 
made before the summer of 2006, and has since been fixed at the 
GMRT. Initial editing of the data was performed separately on each 
sideband with standard AIPS tasks, to remove bad baselines, anten- 
nas, and channels that were suffering from large amounts of narrow 
band interference, along with the first and last integration periods 
of each scan. The flux calibrators were used to create a bandpass 
correction for each antenna. In order to create a continuum chan- 
nel, five central frequency channels were combined together, and an 
antenna-based phase and amplitude calibration created using obser- 
vations of J1634+627. This calibration was applied to the original 
data, which was then compressed into 1 1 channels, each with band- 
width of 1.25 MHz (so the first few and last few spectral channels 
were omitted from the data, since they tended to be the noisiest). 
The new channel width is small enough that bandwidth smearing 
is not a problem in our images, and led to an effective bandwidth 
of 13.75 MHz in each sideband. Further flagging was performed 
on the 11 channel dat a set, and the tw o sidebands combined us- 
ing UVFLP (again, see iGarn et al.ll2007 ^ to improve the uv cover- 
age. The coverage for one of the shallow observations is shown in 
Fig.|2] Baselines shorter than 1 kA were omitted from the imag- 
ing, since the GMRT has a large number of small baselines which 
would otherwise dominate the beam shape. 

Each pointing was broken down into 3 1 smaller facets (as dis- 
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cussed in lGam et al ] [2007h . arranged in a hexagonal grid and cov- 
ering an area with diameter ~ 1?8. These were imaged separately, 
each with a different assumed phase centre. The large area cov- 
ered (compared with the Full-Width Half-Maximum of the GMRT, 
which is ~ 0?74) allows bright sources well outside of the observed 
region to be cleaned from the images, while the faceting procedure 
avoids the introduction of phase errors due to the non-planar nature 
of the sky. All images were made with the same elliptical synthe- 
sised beam of size 6x5 arcsec^, position angle +45° by setting the 
parameters for the restoring beam within IMAGR. A pixel size of 
1.5 arcsec was chosen, to ensure that the beam was well oversam- 
pled. 

Each pointing went through three iterations of phase self- 
calibration at 10, 3 and 1 minute intervals, and then a final round 
of self-calibration correcting both phase and amplitude errors, at 
10 minute intervals. The overall amplitude gain was held constant 
in order not to alter the flux density of sources. The self-calibration 
steps improved the noise level by about 10 per cent, and signifi- 
cantly reduced the residual sidelobes around the brighter sources. 

The four deep pointings have a final rms of between 40 and 
42 /iJy before correction for the GMRT primary beam, while the 
15 shallow fields have a noise level of between 66 and 73 /xJy. 
These figures are very close to the expected thermal noise limits of 
36 and 63 /jJy respectively, calculated using 



G\Jn{n — 1)A'ifAvt 



(1) 



where the system temperature ~ 92 K, and the antenna gain G ~ 
0.32 K Jy"' - values taken from the GMRT websittQ — « is the 
number of working antennas, typically 27 during our observations, 
N\Y = 2 is the number of sidebands, Av = 13.75 MHz is the effec- 
tive bandwidth per sideband, and T is the integration time for each 
pointing; 

In bam et alj ( l2007h . we detected a position-dependent error 
with the GMRT primary beam. This led to a systematic difference 
between the measured flux densities of sources that were present 
in more than one pointing, with the fractional offset varying across 
the primary beam. To check for a similar effect here, we corrected 
each pointing for the nominal beam s hape, using an 8th - order poly- 
nomial with coefficients taken from lKantharia & Raol ( 1200 ih . We 
took the 557 sources present in the overlapping regions between 
two pointings with peak brightness greater than 1.5 mjy beam^', 
and looked at the fractional offsets as a function of position away 
from each pointing centre. Fig. [3] shows the variation in this offset 
with source Right Ascension. A similar trend was seen with Decli- 
nation. 

We repeated the analysis of iGam et 'Zl j2007l) , to model the 
effective pointing centre of the telescope as having a systematic 
offset from its nominal value. By shifting the phase centre of all 
pointings by ~ 2f7 in a north-west direction before performing the 
primary beam correction, we were able to remove this systematic 
effect, as shown in Fig|4l The amount and direction of the correc- 
tion was consistent between pointings, and did not vary between 
the deep and shallow observations. This correction also removed 
the systematic effect seen in the Declination direction. The size and 
direction of this correction is similar to that seen in the xFLS sur- 
vey. 

The accuracy of the primary beam correction was then tested, 
as shown in Fig. |5] The fractional offset of sources in overlapping 
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Figure 3. Fractional offset in peak brightness of all sources in overlap- 
ping regions, using the nominal primary beam centre. A gradient across 
the pointing is seen, leading to systematic errors in measured flux density. 
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Figure 4. Fractional offset in peak brightness of sources in overlapping 
regions, using the shifted primary beam centre applied to all fields. 



regions is now a good fit to a Gaussian, with mean of 0.009 and 
a of 0.1, indicating a 10 per cent error in the absolute flux density 
calibration of all sources. 

The 19 pointings were mosaicked together, taking into ac- 
count the offset primary beam and weighting the final mosaic ap- 
propriately by the relative noise of each pointing. The mosaic was 
cut off at the point where the primary beam correction dropped to 
20 per cent of its central value. The final mosaicked image is avail- 
able via http : / /www . mrao . cam . ac. uk/ surveys/ 
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Figure 5. Distribution of fractional offsets, after correction for the shifted 
primary beam. The best-fit Gaussian is shown, giving a 10 per cent error in 
the overall flux density calibration. 
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Figure 6. The rms noise of the final mosaic, calculated using Source Ex- 
tractor. The grey-scale ranges between 40 and 350 iliy, and the contours 
are at 60 and 120 ;tlJy respectively. 

3 RESULTS AND DISCUSSION 

Source Extractor ("SExtractor jBertin & Arnoutill996h was used to 
calculate the rms noise CT across the mosaic. A grid of 16 x 16 pixels 
was used in order to track changes in the local noise level, which 
varies significantly near the brightest sources. Fig.[6]illustrates the 
local noise, with the grey-scale varying between 40 ;tiJy (the noise 
level in the centre of the deep pointings) and 350 /ijy (the noise 
level for the shallow pointings, at the distance where the GMRT 
primary beam gain was 20 per cent of its central value). The 60 and 
120 /im contours are plotted, which cover the majority of the deep 
and shallow survey regions respectively. 

A sample region of the ELAIS-Nl survey is shown in Fig.|7]to 
illustrate the quality of the image. Most of the sources in our survey 
are unresolved, although there are several with extended structures. 
We present a sample of these in Fig. [8] 

Our GMRT data suffer from dynamic range problems near 
the brightest sources, and the final mosaic has increased noise and 
residual sidelobes in these regions. We had fewer problems with 
our survey of the xFLS region, due to the longer time spent on 
each pointing and the correspondingly better uv coverage. There 
were also fewer bright sources in the xFLS field, so a much smaller 
region was affected by residual sidelobes. Fig. [9] shows an area 
around one of the bright sources, to illustrate the problems caused 
by the residual sidelobes. While the local noise calculated by SEx- 
tractor increases due to these residuals, some of them still have an 
apparent signal-to-noise level that is greater than 6. We therefore 
opted for a two-stage selection criteria for our final catalogue. 

3.1 Source fitting 

An initial catalogue of 4767 sources was created using SExtrac- 
tor. The mosaic was cut off at the point where the primary beam 
correction dropped to 20 per cent of its central value (a radius of 
0.°53 from the outer pointings), however only sources inside the 
30 per cent region (0?47) were included in the catalogue to avoid 
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Figure 9. Source GMRTENl J161212.4+552308, and the errors suiTound- 
ing it. The grey-scale ranges between —0.2 and 1 mJy beam^', and the 
source has a peak brightness of 389 mJy. The region affected by an over- 
density of sources is shown by the black circle, of radius 10 arcmin - see 
text for more details. 

the mosaic edges from affecting the estimation of local noise. The 
requirements for a source to be included were that it had at least 
5 connected pixels with brightness greater than 3c7, and a peak 
brightness greater than 6cT. The image pixel size meant that the 
beam was reasonably oversampled, so the source peak was taken 
to be the value of the brightest pixel within a source. The inte- 
grated flux density was calculated using the FLUX_AUTO option 
within SExtractor. This creates an e lliptical aperture around each 
object (as described in lKronlll980h . and integrates the flux con- 
tained within the ellipse. Comparisons between the flux density 
obtained through this method and through the method developed in 
iGam et al. I j2007h - pixels above a given threshold were summed, 
then empirically corrected for the elliptical beam shape - give good 
agreement between the two techniques. Sources with Kron flux 
density above 1 mJy showed no statistical difference between the 
two flux density measurements, with an uncertainty of 3 per cent. 
Sources below 1 mJ y had a Kron flux density that was systemati- 
cally larger than the lGam et al.l ( |2007|) method by 6 per cent, with 
an uncertainty of 4 per cent. We chose to use the Kron method in 
order to avoid the empirical correction factor. 

In order to estimate the area affected by artefacts near the 
bright sources, we calculated the number of (potentially spurious) 
sources in a series of concentric rings centred on a bright source, 
and converted this to an effective source density as a function of 
distance. While the noise does vary across the map, the variation is 
smooth and relatively slow away from the bright sources. Using the 
more distant rings we calculated a mean source sky density /i, and 
an estimate of the error in this value, Gfi, which will be unaffected 
by the presence of the spurious sources, and will vary slightly be- 
tween each bright central source due to the changing properties of 
the image. 

Fig.[lO]shows the source density (in arbitrary units) away from 
the source in Fig.|9] The clear over-density of sources can be seen 
within 10 arcmin of the source. If there is a significant peak in the 
density (greater than /I + 6ct^), then we define the size of the af- 
fected region by finding the first radius at which the source density 
drops below /x + 3(7^1 - for this source, the affected radius is 10 ar- 
cmin. The affected region has been added to Fig. [9] Within this 
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Figure 7. A 70 x 70 arcmin^ region of the 610 MHz image, to illustrate the quality of the survey. The region is located within the deeper area of the survey, 
and most sources are unresolved. The grey-scale ranges between —0.2 and 1 mJy beam^', and the noise is relatively uniform and between 40 and 60 ^Jy, 
apart from small regions near bright sources where the noise increases. 



region, only sources with a peak brightness greater than 12c7 are 
included in our catalogue - this value was determined empirically. 
The source density plot for a 10 mJy source is shown in Fig.[TT] on 
the same scale as Fig. 1 101- while there may still be a slight overden- 
sity near the centre, the increased noise level in this region, along 
with the 6(7 cut-off reduces the risk of selecting a spurious source 
near sources with weak over-densities. 

This analysis was repeated for all sources with a peak greater 
than 10 mJy in order to filter spurious sources. The final catalogue 
contains 2500 sources - we have erred on the side of caution in or- 
der to produce a catalogue with little contamination from spurious 
sources. The size of the affected region is correlated with the peak 
brightness of a source (with Pearson product-moment correlation 
coefficient of 0.53), and the number of spurious sources is also cor- 
related with the peak brightness, with correlation coefficient 0.73. 
The precise size of the affected region depends on the uv coverage 



for the relevant pointing, the time spent on observations and the 
local noise levels. 

Table 13.11 presents a sample of 60 entries in the cata- 
logue, which is sorted by right ascension. The full table is avail- 
able via |http : / /www . mrao . cam . ac. uk/ surveys/ 1 Column 1 
gives the lAtj designation of the source, in the form GM- 
RTENl Jhhmmss.s+ddmmss, where J represents J2000.0 coordi- 
nates, hhmmss.s represents right ascension in hours, minutes and 
truncated tenths of seconds, and ddmmss represents the declina- 
tion in degrees, arcminutes and truncated arcseconds. Columns 2 
and 3 give the right ascension and declination of the source, calcu- 
lated by first moments of the relevant pixel brightnesses to give a 
centroid position. Column 4 gives the brightness of the peak pixel 
in each source, in mJy beam^', and column 5 gives the local rms 
noise in /ijy beam"'. Columns 6 and 7 give the integrated flux 
density and error, calculated from the local noise level and source 
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Figure 8. A selection of extended objects in the ELAIS-Nl 610 MHz GMRT survey - contours are plotted at ±200 ^iJy x 1, V2, 2, 2V2, 4. . ., with the 
exception of object (e), where the contours start at ±500 ^Jy. Negative contours are represented by dashed lines. The resolution of the beam is shown in the 
bottom left of each image, and the designations of each source component are given below. 



A 610-MHz survey of the ELAIS-Nl field 7 



s 

Q 



0.004 



0.002 




20 

Radius (arcmin) 



40 



Figure 10. Density of sources, in concentric rings around the bright source 
shown in Fig. [9] The overdensity of sources near to the bright (389 mjy) 
central object can be clearly seen. The mean source density, far from the 
central source, is given by the large dashed line while the cutoff density 
defining the affected region is given by short dashes - see text for more 
details. 
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Figure 11. Density of sources, in concentric rings around a fainter (10 mJy) 
object in the ELAIS-Nl survey field. There is still a slight overdensity near 
the central source, but the signal-to-noise cutoff of 6(T ensures that spurious 
sources are not included in the catalogue. 



size. Columns 8 and 9 give the X, Y pixel coordinates of the source 
centroid from the mosaic image. Column 10 is the Source Extrac- 
tor deblended object flag - 1 where a nearby bright source may be 
affecting the calculated flux, 2 where a source has been deblended 
into two or more components from a single initial island of flux, 
and 3 when both of the above criteria apply. There are 232 sources 
present in our catalogue with non-zero deblend flags; it is necessary 
to examine the images to distinguish between the case where one 
extended object has been represented by two or more entries, and 
where two astronomically distinct objects are present. 



3.2 Comparison with other surveys 

In order to test the positional accuracy of our catalogue, we paired 
it with the 393 obj ects in the VLA 1.4 GHz source catalogue of 
ICiliegietal.Ul999ll) . using a pairing radius of 6 arcsec. The VLA 
survey covers only ~ 25 per cent of our 610 MHz observations. 
Fig.ll2l shows the position offset of the 263 matched sources com- 
pared with their VLA counterparts - the distribution of offsets is 
approximately Gaussian, with mean offset in Right Ascension of 
—0.9 arcsec, and —0. 1 arcsec in Declination. The standard devia- 
tions of the distribution are 0.9 and 0.7 arcsec respectively. 

We repeated th is analysis with data from the FIRST survey 
( lBeckeretai]|l995[) . The whole of our 610 MHz survey region 
is covered by FIRST, and even with the reduced sensitivity of 
FIRST (--^150 /xJy noise), 504 pairs are found within 6 arcsec. 
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Figure 12. Source positions in th e GMRT catalogue r elative to the positions 
found in the VLA catalogue of ICiliegi et alj i 19991) . for unique matches 
within 6 arcsec. Offsets are distributed in a Gaussian fashion in RA and 
DEC, and the ellipse corresponding to 1 fj for the distribution is shown. 
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Figure 13. Radio spectral index a between 610 MHz and 1.4 GHz, for 
sources in the VLA ELAIS-Nl catalogue of Ciliegi et al. (solid lines) and 
in the FIRST catalogue (dashed lines). 



We again find Gaussian-distributed position errors, with RA off- 
set of —0. 1 arcsec, standard deviation 0.4 arcsec and DEC offset of 
—0. 1 arcsec, standard deviation 0.6 arcsec. We have not corrected 
the positions given in our catalogue, since it agrees closely with 
FIRST. 

The spectral index distribution of the matched sources from 
both surveys is shown in Fig. [13] using the integrated flux density 
measurements. The distribution peaks around a = 0.7, where a 
is defined such that the flux density 5 scales with frequency v as 
5 = 5oV-«. 

Fig. [14] shows the spectral index distribution for all sources 
with matches in the GMRT and FIRST catalogues (black diagonal 
crosses), and sources in the GMRT and Ciliegi et al. ( 1999) cata- 
logues (red upright crosses). There are significant biases in Fig. 1141 
due to the varyin g sensitivity l evels o f the three surveys. In the re- 
gion covered by ICiliegi et 

m jl999l) , the 610 MHz completeness 
level is 360 ^Jy, shown by the black dotted line. The limiting spec- 
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Table 1. A sample of 60 entries from the 610-MHz ELAIS-Nl catalogue, sorted by right ascension. The full version of this table is available as Supplementary 
Material through the online version of this article, and via http : //www . mrao . cam . ac . uk/ surveys/ 
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Int. Flux Density 


Error 


X 


Y 


Flags 






J2000.0 


J2000.0 


mJy beam^ ' 


/Jjy beam"' 


mJy 


mJy 










(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


GMRTENl 


J160319.2+542543 


16:03:19.22 


+54:25:43.6 


0.496 


75 


0.415 


0.067 


6599 


2946 





GMRTENl 


J160319.2+553149 


16:03:19.24 


+55:31:49.1 


0.469 


66 


0.383 


0.054 


6527 


5589 





GMRTENl 


J160319.3+554950 


16:03:19.38 


+55:49:50.2 


0.514 


85 


0.265 


0.054 


6506 


6309 





GMRTENl 


J160320.8+550645 


16:03:20.89 


+55:06:45.9 


2.617 


97 


3.425 


0.175 


6545 


4587 





GMRTENl 


J160321. 1+553654 


16:03:21.13 


+55:36:54.2 


0.418 


69 


0.309 


0.076 


6510 


5792 





GMRTENl 


J160322.6+554320 


16:03:22.67 


+55:43:20.8 


0.560 


75 


0.444 


0.068 


6495 


6049 





GMRTENl 


J 160323. 1+543737 


16:03:23.13 


+54:37:37.3 


1.920 


74 


2.256 


0.119 


6564 


3421 





GMRTENl 


J160324.0+550618 


16:03:24.05 


+55:06:18.2 


0.635 


84 


0.495 


0.079 


6527 


4568 





GMRTENl 


J160326.0+540817 


16:03:26.00 


+54:08:17.1 


0.595 


94 


0.486 


0.081 


6579 


2248 





GMRTENl 


J160327.7+552647 


16:03:27.78 


+55:26:47.3 


3.240 


84 


3.867 


0.149 


6484 


5387 





GMRTENl 


J160327.9+543326 


16:03:27.99 


+54:33:26.0 


0.452 


75 


0.338 


0.064 


6540 


3253 





GMRTENl 


J160329.5+540705 


16:03:29.55 


+54:07:05.7 


0.869 


131 


0.681 


0.117 


6559 


2200 





GMRTENl 


J160330.8+542454 


16:03:30.83 


+54:24:54.3 


0.483 


78 


0.837 


0.097 


6533 


2912 





GMRTENl 


J160331. 1+554327 


16:03:31.14 


+55:43:27.2 


0.473 


76 


0.555 


0.071 


6447 


6052 





GMRTENl 


J160331.3+545000 


16:03:31.39 


+54:50:00.7 


0.408 


67 


0.322 


0.055 


6503 


3916 





GMRTENl 


J160332.3+553000 


16:03:32.38 


+55:30:00.2 


0.508 


69 


0.536 


0.076 


6454 


5514 





GMRTENl 


J160332.6+554622 


16:03:32.66 


+55:46:22.6 


3.212 


76 


4.152 


0.153 


6435 


6169 





GMRTENl 


J160332.9+541746 


16:03:32.90 


+54:17:46.8 


0.474 


73 


0.381 


0.060 


6528 


2627 





GMRTENl 


J160333. 1+542914 


16:03:33.11 


+54:29:14.2 


2.377 


73 


9.697 


0.226 


6515 


3085 


3 


GMRTENl 


J160333.6+552623 


16:03:33.69 


+55:26:23.7 


0.773 


78 


0.878 


0.101 


6451 


5370 





GMRTENl 


J160333.7+540540 


16:03:33.70 


+54:05:40.6 


7.541 


294 


9.799 


0.552 


6536 


2142 





GMRTENl 


J160334.6+542900 


16:03:34.68 


+54:29:00.2 


2.326 


80 


4.235 


0.167 


6506 


3075 


3 


GMRTENl 


J160335. 1+551534 


16:03:35.19 


+55:15:34.1 


0.462 


74 


0.333 


0.066 


6454 


4937 





GMRTENl 


J160335.2+555419 


16:03:35.24 


+55:54:19.5 


0.762 


105 


0.984 


0.108 


6412 


6486 





GMRTENl 


J160335.2+540515 


16:03:35.25 


+54:05:15.3 


38.494 


383 


50.351 


0.961 


6528 


2125 





GMRTENl 


J160336.5+555147 


16:03:36.55 


+55:51:47.5 


0.500 


79 


0.662 


0.082 


6408 


6385 





GMRTENl 


J160336.9+544120 


16:03:36.90 


+54:41:20.7 


0.532 


68 


0.853 


0.092 


6480 


3568 





GMRTENl 


J160337.6+545944 


16:03:37.64 


+54:59:44.2 


0.547 


80 


0.530 


0.077 


6457 


4303 





GMRTENl 


J160338.7+554348 


16:03:38.77 


+55:43:48.1 


0.771 


80 


1.743 


0.135 


6404 


6065 





GMRTENl 


J160339.0+545943 


16:03:39.05 


+54:59:43.8 


0.500 


80 


0.415 


0.065 


6449 


4303 





GMRTENl 


J160339.3+551352 


16:03:39.37 


+55:13:52.0 


0.516 


78 


0.414 


0.063 


6432 


4868 





GMRTENl 


J160339.6+542953 


16:03:39.68 


+54:29:53.2 


0.701 


66 


0.647 


0.078 


6476 


3110 





GMRTENl 


J160339.7+550600 


16:03:39.75 


+55:06:00.4 


0.584 


80 


0.418 


0.069 


6438 


4554 





GMRTENl 


J 160340. 1+545 127 


16:03:40.10 


+54:51:27.4 


0.711 


79 


0.639 


0.079 


6451 


3972 





GMRTENl 


J160340.1+550544 


16:03:40.15 


+55:05:44.9 


0.495 


81 


0.392 


0.069 


6436 


4543 





GMRTENl 


J160340.8+554325 


16:03:40.83 


+55:43:25.8 


1.447 


81 


4.633 


0.210 


6392 


6050 





GMRTENl 


J160341. 2+552611 


16:03:41.27 


+55:26:11.8 


0.475 


68 


0.460 


0.065 


6408 


5361 





GMRTENl 


J160341.5+552205 


16:03:41.59 


+55:22:05.4 


0.638 


69 


0.614 


0.080 


6411 


5197 





GMRTENl 


J160341.6+553203 


16:03:41.66 


+55:32:03.4 


0.462 


69 


0.365 


0.062 


6400 


5595 





GMRTENl 


J 160343. 1+540324 


16:03:43.13 


+54:03:24.6 


0.908 


126 


0.556 


0.092 


6483 


2050 





GMRTENl 


J 160344. 8+540320 


16:03:44.84 


+54:03:20.4 


1.168 


118 


0.856 


0.110 


6473 


2047 





GMRTENl 


J160345. 8+553021 


16:03:45.84 


+55:30:21.4 


0.381 


60 


0.370 


0.065 


6378 


5527 





GMRTENl 


J160345. 8+554238 


16:03:45.88 


+55:42:38.3 


0.653 


76 


2.566 


0.156 


6365 


6018 





GMRTENl 


J160346.5+552855 


16:03:46.56 


+55:28:55.8 


0.577 


73 


0.510 


0.076 


6375 


5469 





GMRTENl 


J160346.6+550826 


16:03:46.62 


+55:08:26.0 


0.525 


70 


0.416 


0.065 


6396 


4650 





GMRTENl 


J160348.3+542626 


16:03:48.37 


+54:26:26.4 


2.328 


80 


2.660 


0.133 


6429 


2970 





GMRTENl 


J160348.6+550124 


16:03:48.61 


+55:01:24.3 


0.477 


78 


0.257 


0.050 


6392 


4368 





GMRTENl 


J160349.2+554243 


16:03:49.21 


+55:42:43.9 


3.349 


85 


4.224 


0.163 


6346 


6021 





GMRTENl 


J160350.0+545634 


16:03:50.06 


+54:56:34.6 


0.435 


67 


0.371 


0.058 


6389 


4175 





GMRTENl 


J160350.4+550717 


16:03:50.43 


+55:07:17.2 


0.597 


78 


0.626 


0.083 


6376 


4603 





GMRTENl 


J160350.5+541302 


16:03:50.51 


+54:13:02.7 


0.770 


80 


0.535 


0.068 


6430 


2435 





GMRTENl 


J160351. 1+555644 


16:03:51.13 


+55:56:44.4 


3.996 


120 


7.391 


0.287 


6321 


6581 





GMRTENl 


J16035 1.2+552346 


16:03:51.24 


+55:23:46.1 


0.519 


80 


0.640 


0.072 


6354 


5262 





GMRTENl 


J16035 1.3+540609 


16:03:51.38 


+54:06:09.9 


0.951 


110 


0.843 


0.110 


6432 


2159 





GMRTENl 


J160352.6+550012 


16:03:52.62 


+55:00:12.8 


0.510 


71 


0.592 


0.076 


6370 


4320 





GMRTENl 


J160352.8+542943 


16:03:52.80 


+54:29:43.5 


2.258 


82 


2.422 


0.124 


6400 


3101 





GMRTENl 


J160353.2+550309 


16:03:53.29 


+55:03:09.9 


0.570 


74 


0.500 


0.072 


6363 


4438 





GMRTENl 


J160355.5+553844 


16:03:55.52 


+55:38:44.9 


0.463 


73 


0.242 


0.050 


6314 


5861 





GMRTENl 


J160356.2+550126 


16:03:56.20 


+55:01:26.7 


0.472 


77 


0.380 


0.063 


6348 


4369 





GMRTENl 


J160356.2+550415 


16:03:56.25 


+55:04:15.3 


0.506 


83 


0.451 


0.068 


6345 


4481 
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Figure 14. The variation in spectral index a with 610 MHz flux den- 
sity. Black diagonal crosses represent sources in the GMRT and FIRST 
catalogues, with the solid black line showing the limiting spectral index 
that could be detected, given the respective sensitivity levels. Red upright 
crosses represent sources in the corresponding GMRT and Ciliegi et al. cat- 
alogues, with the dashed red line showing the limit on a. The 610 MHz flux 
density hmit is shown by the dotted black line. 



tral indices for sources at the sensitivity levels of the 1.4 GHz sur- 
veys are shown (FIRST ~ solid black line, Ciliegi - dashed red line). 
In order to look for variations in the source population, we calcu- 
late the mean and median spectral indices for sources with detec- 
tions in the Ciliegi catalogue, with 610 MHz flux density between 
500 /iJy and 1 mJy - the point at which the turnover in source 
counts becomes visible (see Section |4l( ~ and above 1 mJy. The 
mean values of a are 0.22 ± 0.09 and 0.45 ± 0.04 respectively, and 
the median values are 0.36 ±0.12 and 0.56 ±0.04. There are 48 
and 168 sources in the two flux density bins. The bias against steep- 
spectrum sources at low flux densities (which is visible in Fig.ll4t 
means that, for a source with 610 MHz flux density of 500 /iJy, 
the largest value of a that would be detectable is -^0.8 and so this 
apparent flattening at f ainter flux densities may simply be due to 
sample bias. However, iBondi et alj j2007h also find significantly 
flatter spectral indices for fainter radio sources, again comparing 
610 MHz and 1.4 GHz data, and attribute this to the emergence 
of a population of low-luminosity AGNs - see Section|4]for more 
details. 



4 DIFFERENTIAL SOURCE COUNTS 

We derived source counts for the ELAIS-Nl field by considering 
two separate regions of the map. In order to study the 610 MHz 
source counts at low flux densities, we selected a 70 x 70 arcmin^ 
region (shown in Fig. |7}, which was free from spurious sources 
above the local 6C7 value. The region comes from the deeper sec- 
tion of the survey field. Fig.[T5]shows the distribution of pixel noise 
levels (taken from the SExtractor rms map), along with the cumu- 
lative fraction of pixels with noise below a given flux density. 

We constructed source counts by binning our sources by their 
integrated flux density, with bins ranging between 270 fxJy and 
10 mJy. The lower limit was selected as the 6<J flux density cor- 




60 70 
Noise (microJy) 

Figure 15. The fraction of pixels in the 70 x 70 arcmin^ region of the SEx- 
tractor rms map with each noise level (boxes, and left-hand axis), along with 
the cumulative fraction of pixels with noise below a particular value (sohd 
line, and right-hand axis). 



responding to a noise of 45 /xJy, close to the peak in noise dis- 
tribution seen in Fig. [15] We selected a relatively low upper limit 
of 10 mJy, since our region was specifically chosen for its lack of 
brighter sources (and the errors resulting from these sources). 

The source counts were corrected for the fraction of the image 
over which they could be detected, taking into account the increase 
in noise near the bright sources, and for the resolution bias inherent 
in constructing flux density counts from a catalogue with a signal- 
to-noise ratio based around the peak brightness of a source. We 
use the resolution corrections of iMoss et al. I|2007), who found that 
~7 per cent of sources below 360 fiJy and 3 per cent of sources 
between 360 and 1000 /xJy were missed from a similar resolu- 
tion GMRT 610 MHz survey of the XMM-Newton/Chandra sur- 
vey field. No cor rection for Eddington bias has been made, since 
IMoss et all ( 120071) find that it only affects their faintest flux density 
bin, increasing it by approximately 20 per cent. 

The differential source count dA' / d5 was calculated by divid- 
ing Nc, the corrected number of sources in each bin, by AA5, where 
A is the total area of the image in steradians, and hS is the width of 
the flux bin in Jy, i.e. 



(2) 



dA' _ A'c 
dS " AA5 

We calculated the differential source count for sources with 
flux densities above 10 mJy by considering the central --^ 4 deg^ of 
the map. This region has several bright sources with residual side- 
lobes, but inspection of these residuals found them to have typical 
flux densities of below 1 mJy, which will therefore not affect the 
bright source counts. No corrections for resolution bias or Edding- 
ton bias have been made for these brighter sources. 

Table |4] gives the source counts, mean flux density (5) of 
sources in each bin, d/V/d5 and d/V/d5 normalised by (5)^'^, the 
value expected from a static Euclidean universe. dA'/d5 is plot- 
ted in Fig. [161 along with 610 MHz source counts from fur- 
ther GMRT surv eys of the XMM-Newton/Ch andra survey fiel d 
jMoss et al.120071) , the VVDS-VLA Deep Field dSondi et al.l2007h . 
and our previous survey of the xFLS field (Gam et al. l2007h . A se- 
ries of Wester bork Synthesis Radio Telesc ope (W SRT) surveys at 
610 MHz (bylValentiin et aL|[l977l ; lKatgerti[l979l : IValentiin|[T980l; 
iKatgert-MerkeUin et al.lll985l) are also shown for comparison. 

Fig.lI7lshows the normalised differential source counts. There 
is good agreement in the various surveys above 1 mJy, and all sur- 
veys show a turnover just above 1 mJy although there is approx- 
imately a factor of 2 discrepancy between the highest and lowest 
points in the value of dN/dS (S)^'^ below this flux density. A varia- 
tion in_source_countsJias been previously observed at 1.4 GHz (se e 
e.g. lSevmour et alj2004lHuvnh et al.l2005l : lBiggs & Ivisonl2006h . 
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Figure 16. Differential source counts 4^ at 610 MHz from surveys taken with the WSRT iValentiin et all 1 19771 ; iKatgerll Il979t IValentiinI Il980l ; 
[Katgert-Merkeliin et al...l985.^ and GMRT ffiondi et al. 2007 : Moss et al. 2007; Garn et al. 200?!) 



Table 2. 610 MHz differential source counts for the ELAIS-Nl survey. 



Flux Bin 
(raJy) 



(S) 
(raJy) 



N 



dN/dS 
(sr-l Jy-l) 



dN/dS {S)--^ 
(sr-' Jy'-5) 



0.270 - 


0.387 


0.380 


89 


195.9 


4.0±0.4x 


10' 


8.44±0.89 


0.387 ^ 


0.556 


0.462 


100 


111.9 


1.6±0.2x 


lO' 


7.34±0.73 


0.556 - 


0.798 


0.665 


67 


69.2 


6.9±0.8x 


10^ 


7.87 ±0.96 


0.798 - 


1.145 


0.953 


45 


45.0 


3.1±0.5x 


10^ 


8.77 ± 1.31 


1.145- 


1.643 


1.396 


43 


43.0 


2.1±0.3x 


10« 


15.2±2.31 


1.643- 


2.358 


1.950 


24 


24.0 


8.1±1.6x 


10^ 


13.6±2.77 


2.358- 


3.384 


2.775 


24 


24.0 


5.6±1.2x 


10^ 


22.9 ±4.67 


3.384- 


4.856 


3.989 


15 


15.0 


2.4±0.6x 


10^ 


24.7 ±6.37 


4.856- 


6.968 


5.971 


13 


13.0 


1.5±0.4x 


10^ 


40.9 ±11.3 


6.968 - 


10.00 


8.394 


13 


13.0 


1.0±0.3x 


10^ 


66.7±18.5 


10.00 - 


13.49 


11.41 


33 


22.0 


8.1±1.4x 


10*^ 


112±19.6 


13.49- 


18.20 


15.73 


19 


19.0 


3.4±0.8x 


10^ 


107 ±24.6 


18.20- 


24.56 


21.93 


18 


18.0 


2.4±0.6x 


10* 


173 ±40.6 


24.56 - 


33.14 


29.41 


14 


14.0 


1.4±0.4x 


10* 


207 ±55.3 


33.14- 


44.72 


38.92 


10 


10.0 


7.4±0.2x 


105 


221 ±69.8 


44.72 - 


60.34 


53.74 


2 


2.0 


l.liO.l X 


10^ 


73.3±51.9 


60.34 - 


81.41 


68.60 


3 


3.0 


1.2±0.7x 


105 


150.1 ±86.7 


81.41 - 


109.8 


00.00 





0.0 


0.0 




0.0 


109.8 - 


148.2 


134.0 


1 


1.0 


2.2 ± 2.2 X 


10-* 


147 ±147 


148.2- 


200.0 


191.9 


1 


1.0 


1.7±1.7x 


10" 


267 ±267 



In particular, iBiggs & Ivisonl ( l2006h and lSimpson et alj l l2006l) find 
a factor of 2 difference in the 1 .4 GHz source counts for different 
fields at ~180 /iJy, equivalent to ~350 /xJy at 610 MHz for an as- 
sumed spectral index of 0.8. This difference could be due to cosmic 
variance and the small numbers of sources detected in each field at 
these depths, and further data is needed to determine where the true 
source counts lie at these flux density levels. 

The turnover in source counts is a well-known feature of the 
source counts at 1.4 GHz, and is thought to be due to the emer- 
gence o f a new po pulation of sources, with some authors (e.g . 
Condon|[l9 89; Ro wan-Robinson et al.|[l993l ; iHopkins et alj|l99i : 
Seymour et al.ll2004l) attributing it to star-forming galaxies, while 



others jjarvis & Rawlingsl2004l : [simDson et al.l2006l) claiming that 
low-luminosity AGNs may form an important contribution to the 
source counts at these levels. The turnover occurs at 1 mjy at 
1.4 GHz, equivalent to --^1.9 mJy at 610 MHz for a spectral index 
of 0.8, which is consistent with Fig. 1171 

In Fig. [18] we show the differential source counts from the 
combination of our ELAIS-Nl and Spitzer extragalactic First Look 
Survey field - there are more than a hundred sources in the low- 
est four flux bins, leading to the most precise estimation of the 
610 MHz source counts to date. 
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Figure 18. The normalised 610 MHz differential source counts, found by combining this work and the xFLS survey of lGarn et alj jlOOTi) . 
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